Discussion on exp-function method and 
modified method of simplest equation 



Zlatinka I. Dimitrova 



Abstract 

We discuss the relation between the modified method of simplest 
equation and the exp-function method. First on the basis of our ex- 
perience from the application of the method of simplest equation we 
generalize the exp-function ansatz. Then we apply the ansatz for 
obtaining exact solutions for members of a class of nonlinear PDEs 
which contains as particular cases several nonlinear PDEs that model 
the propagation of water waves. 
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1 Nonlinear PDEs and method of simplest 
equation 

Nonlinear models are much used in various branches of science f 1 ' 2,3 ]. Often 
such models contain nonlinear PDEs and because of this the interest in ob- 
taining exact analytical solutions of nonlinear PDEs increses steadily. Such 
exact solutions often describe important classes of waves and processes in 
the investigated systems. In addition the exact solutions can be useful as 
initial conditions in the process of obtaining of numerical solutions or as test 
soluitions for computer programs for obtaining numerical solutions of the 
studied nonlinear PDEs. 

Because of all above the nonlinear PDEs are widely applied in the theory 
of solitons [ 4,s ], hydrodynamics and theory of turbulence [ 6 ] - [ 10 ], theory of 
dynamical systems, chaos [ u ] - [ 14 ], etc. Sophisticated methods for obtaining 
exact solutions of nonlinear PDEs such as the inverse scattering transform or 
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the method of Hirota [ 15 ] allow obtaining of soliton solutions of some equa- 
tions. In the last several years effective approaches for obtaining exact special 
solutions of complicated nonlinear nonintegrable PDEs have been developed 
too [ 16 ] - [ 18 ]. These approaches leaded to exact solutions of many equations 
such as the Kuramoto-Shivasinsky equation [ 19 ] or equations, connected to 
the models of migration of populations[ 20 ' 21 ' 22 ]. The discussion below will be 
devoted to the modified method of simplest equation: a version of the method 
of simplest equation for obtaining exact solutions of nonlinear PDEs and on 
the relation of this method to another popular method: the exp-function 
method. 

A brief description of the method of simplest equation is as followsf 23 ]. 
Let us have a partial differential equation and let by means of an appro- 
priate ansatz this equation be reduced to the nonlinear ordinary differential 
equation 



For large class of equations from the kind (pQ) exact solution can be con- 
structed as finite series 



where v > 0, fi, are parameters and is a solution of some ordi- 

nary differential equation referred to as the simplest equation. The simplest 
equation is of lower order than (JTJ) and we know the general solution of the 
simplest equation or we know at least exact analytical particular solution(s) 
of the simplest equation. 

The application of the modified method of simplest equation is as follows. 
First by means of an appropriate ansatz (for an example the traveling-wave 
ansatz) the solved class of nonlinear PDEs is reduced to a class of nonlin- 
ear ODEs of the kind ([[]). In the method of simplest equation the resulting 
ODEs are treated as in the first step of the test for Painleve property: the 
corresponding equation is subject of leading order analysis that leads to de- 
termination of v from Eq.(T5]). In the modified method of simplest equation 
one uses the equivalent procedure of obtaining and solving a balance equa- 
tion as follows. First the finite-series solution ([2]) is substituted in ([T|) and 
as a result a polynomial of $(£) is obtained. Eq. ([2]) is a solution of ([T]) if 
all coefficients of the obtained polynomial of are equal to 0. Then by 
means of a balance equation one ensures that there are at least two terms 
in the coefficient of the largest power of The balance equation gives 
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a relationship between the parameters of the solved class of equations and 
the parameters of the solution. The application of the balance equation and 
setting the coefficients of the polynomial of $(£) to leads to a system of 
nonlinear relationships among the parameters of the solution and the param- 
eters of the solved class of equations. Each solution of the obtained system 
of nonlinear algebraic equations leads to a solution of a nonlinear PDE from 
the investigated class of nonlinear PDEs. 

2 The exp-function method: one possible gen- 
eralization and application 

Let us now consider the exp-function method. The standard exp-function 
ansatz for a solution of a nonlinear partial differential equation is [ 24 ] 

(3) u(x,t) = S = °, - — £ = kx + wt + 5. 

The authors of the exp-function method do not define the class of equa- 
tions for which an exact solution can be obtained by this method. In our 
opinion the method can be applied to some nonintegrable partial differential 
equations with polynomial nonlinearity. The ansatz can be generalized on 
the basis of the following observation. In one of the variants of the modified 
method of simplest equation the one-wave solution of the studied nonlinear 
partial differential equation is searched by the ansatz 



(4) <o = J2mf(0] 



1=0 



where F(£) is a solution of the simplest equation, A\ and B\ are parameters 
(if Bi = I Eq. (j3J) is a polynomial ansatz but B[ can be non-integer number 
too), and £ = x — v t + £ where v is the velocity of the wave and £ is a 
parameter. When the equation of Bernoulli 

dF 

(5) — = aF{0 + b[F(Of\ M = 2,3,..., 

is used as a simplest equation its solutions are 

m = { "^i"- 1 ^} .case a > 0,b < 0, 

[ 1 - 6exp[a(M - + £„)] J 

, , { aexp\a(M - + £„Y| )^ 

6 F(0 = \- T —r^ TWX^i ,casea<0,b>0. 

t l + 6exp[a(M- l)(£ + £o)] J 
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We observe that the term in {. . . } from Eqs. can be easily obtained from 
Eq. ([3]) when n — m—1. 

On the basis of all above the following simple generalization of the exp- 
function ansatz is obtained: One searches for exact solution of the studied 
nonlinear PDE on the basis of the ansatz 



Particular cases of Eq. ([8]) are for an example: (I) The Camassa - Holm 
equation f|-^^+3w§^-2§^|^ — uj^g = 0, that describes the propagation 
of shallow water waves over a flat bottom; (II) The Degasperis - Processi 
equation § - -gfa + Au^ - 3§j*0 - w0 = that is also connected to 
the dynamics of the nonlinear shallow water waves; (III) The Fornberg - 
Whitham equation § - + jjs + ug* - 3^f# - uff = that is used 

^ at ox A ot ox ox ox ox z ox 

as a model for investigation of the wave - breaking. We shall search for 
solutions of the kind (J7|) with Bi = IB where B is a parameter. We shall 
discuss the simplest possible case n = m = 1 . For this case first we shall write 
a balance equation for the maximum powers in the numerators of the terms 
from Eq. ([7]). There are 5 terms in Eq. (j7j). Each term has several powers 
of exp(£) in its numerator. From these several powers one is the maximum 
power for the corresponding term of Eq. ([7]). As the terms in Eq. jTJ) are 
5 we have 5 maximum powers. We impose the balance equation constraint: 
The two largest of these 5 maximum powers must be equal (otherwise some 
parameters of the solution and eventually some parameters of the equation 
will be which in the most of the cases is undesirable) . The equality of the 
two largest powers leads to the balance equation. 

The balance equation constraint leads to the following two possibilities 
for balance equations 




£ = kx + wt + 5 



Let us now apply the above ansatz to the equation 




(9) 




and 



(10) 




4 



We note that L and H are integers and because of this the parameter B 
must have appropriate (even non-integer) values. Let us illustrate this point 
further by means of two small tables for the two possibilities for balance 
equations. 



L 


B 


2LB + 3 


H<2 + ^ 


1 


1 


5 


H < 4 


2 


1 


7 


H < 3 


1 


1/2 


4 


H < 6 


2 


1/2 


5 


H < 4 











Table 1: Several of possible values of the parameters for the case when Eq. 
(jnD is balance equation. 



L 


B 


H = 2 + ^b 


1 


1 


4 


2 


1 


3 


1 


1/2 


6 


2 


1/2 


4 


3 


1/3 


4 









Table 2: Several of possible values of the parameters for the case when Eq. 
(JTUj) is balance equation. 



Let us now consider one examples. Let Eq. fTTOT) be the balance equation 
and in addition B = 1 and L = 1. Then H = 4. Thus the equation we shall 
solve is 

<9w d 3 u d (■sr^ h\ dud 2 u d 3 u 
dt dx 2 dt dx \~^ ahU J 1 ' dxdx 2 ^ U dx 3 
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The ansatz (j7|) becomes 



i 



(12) u(x,t) = ^2A 



i 

1=0 



J2T=0 a i ex P(^) 



ELo 6 i ex p(iO 



, £ = kx + wt + 5 



The substitution of Eq. ( fl2|) in Eq. (iTTjl leads to a system of 5 nonlinear alge- 
braic relationships among the parameters of the solution and the parameters 
of the equation. One solution of this system is 

• a > 0, b < 0: ao = 0; a± = a; bo = 1, b\ — —b 

• a < 0, b > 0: ao = 0; ai = —a; &o = 1; °i = 6 
and 

10A 2 a 2 a 4 - 36ek4o&a 4 ^i - 9A x aba 3 + 18b 2 a 3 A + Qb 2 a 2 + 3Qb 2 a 4 A 2 c 
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k 2 a 2 b 2 

AA\a 2 a 4 - \2aA Q ba 4 Ai - 3A 4 aba 3 + 6b 2 a 3 A + 2b 2 a 2 + Y2b 2 a 4 A\ 



k 2 a 2 b 2 



w = - — — [Aia(a 4 Ala 2 — 8bAiA a 4 a — Aiba 3 a + 6b 2 a 3 A + b 2 a 2 + b 2 a 4 Al) — 
ka l b A 

2A (3b 3 A a 3 - b 3 a 2 - 6b 3 A 2 a 4 )} 

ai = 1 " 2 , 3 [-4A; 2 A\a 4 A ba 4 + 3k 2 A 1 a 3 b 2 a 3 A + 6k 2 l Aia 3 tf l a 4 A 2 + 6b 3 A 2 a 3 + 2b 3 A a 2 + 

l2b 3 A 3 Q a 4 - a A A\a 3 + A\ba 3 a 2 - A 1 ab 2 a 2 + k 2 A\a 5 a 4 + 8bA\A Q a 4 a 2 - QAiab 2 a 3 A - 
18A ia b 2 a 4 A 2 - k 2 A\a%ot 3 + k 2 A x a 3 b 2 a 2 - 3k 2 a 2 A 2 b 3 a 3 - 2k 2 a 2 A b 3 a 2 - Ak 2 a 2 A 3 b 3 a 4 ] 
(13) 

Thus the equation 

du d 3 U du I I r .,;> 1 2 1 , , o,.2 i ,,:i;/2., i C j;l , ,,:!/■> ,2 



+ — [-Ak 2 A 2 a A A ba 4 + 3k 2 Axa 6 b 2 a 3 A + 6k 2 A x a 6 b 2 a 4 A 2 + 



dt dx 2 dt dx \ k 2 a 2 b 3 

6b 3 A 2 a 3 + 2b 3 A a 2 + 12b 3 A 3 a 4 - a 4 A 3 a 3 + A\ba 3 a 2 - A x ab 2 a 2 + k 2 A 3 a 5 a 4 + 
8bA\A a 4 a 2 - 6A 4 ab 2 a 3 A - ^Axatfa^l - k 2 A\a%a 3 + k 2 Aia 3 b 2 a 2 - 

3k 2 a 2 A 2 b 3 a 3 - 2k 2 a 2 A b 3 a 2 - Ak 2 a 2 A 3 Q b 3 a 4 ] + 2a 2 u + 3a 3 u 2 + Aa 4 u 3 J - 

du d 2 u d 3 u 

dx dx 2 ^ dx 3 
(14) 
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has the solutions 



(15) 



<0 



A - det 



A A 1 
a b 



exp(0 



1 - 6exp(0 



a > 0, b < 



A + det 



a b 



exp(0 



(16) 



«(0 



l + 6exp(0 



a < 0, b > 



where 



£ = kx + 



t 



[Ai^aiAla 2 — SbAiAoa^a — Aiba 3 a + 6b 2 a 3 A + b 2 a 2 + tfa^Al) — 

2A (3b 3 A a 3 - b 3 a 2 - 66 3 ^a 4 )] + $ 



ka 2 b 3 



(17) 



and det[. . .] denotes the determinant of the corresponding matrix. 

3 Concluding remarks 

For the nonlinear dynamics, chaos theory and for the nonlinear physics in 
general the methods for obtaining exact analytical solutions of classes of 
nonlinear PDEs are of great interest. It seems however that some of these 
methods are more fundamental than other ones. In this paper we discuss 
the relations between the method of simplest equation and the exp-function 
method. On the basis of our experience gained by application of the method 
of simplest equation to various nonlinear PDEs we consider a generalization 
(171) of the ansatz of the exp-function method and then we demonstrated the 
obtaining of exact traveling wave solutions of a member of the class (JSJ) of 
nonlinear PDEs. 
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